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Abstract — In multiple-antenna communications, as bandwidth 
and modulation order increase, system components must work 
with demanding tolerances. In particular, high resolution and 
high sampling rate analog-to-digital converters (ADCs) are often 
prohibitively challenging to design. Therefore ADCs for such 
applications should be low-resolution. This paper provides new 
insights into the problem of optimal signal detection based on 
quantized received signals for multiple-input multiple-output 
(MIMO) channels. It capitalizes on previous works |1|, |2|, |3|, 
131 which extensively analyzed the unquantized linear vector 
channel using graphical inference methods. In particular, a 
"loopy" belief propagation-like (BP) MIMO detection algorithm, 
operating on quantized data with low complexity, is proposed. 
In addition, we study the impact of finite receiver resolution in 
fading channels in the large-system limit by means of a state 
evolution analysis of the BP algorithm, which refers to the 
limit where the number of transmit and receive antennas go to 
infinity with a fixed ratio. Simulations show that the theoretical 
findings might give accurate results even with moderate number 
of antennas. 



I. Introduction 

Most of the contributions on signal detection for multiple- 
input multiple-output (MIMO) systems assume that the re- 
ceiver has access to the channel data with infinite precision. 
In practice, however, a quantizer (A/D-converter) is applied to 
the received analog signal, so that the channel measurements 
can be processed in the digital domain. In ultra-wideband 
and/or high-speed applications, both the required resolution 
and speed of the ADCs tend to rise, making them expensive, 
power intensive and even infeasible [5|. This work deals with 
MIMO channels with quantized outputs, which we will refer 
to as quantized MIMO systems. 

Recently, graphical inference methods have been applied 
to the usual unquantized MIMO detection problem. In (TJ, 
(2) "loopy" BP-like detection algorithms were derived as low 
complexity heuristics for computing the marginal distribution 
of each signal component. Hereby, we provide an extended 
version of the approximative BP based algorithm relying on 
the commonly used Gaussian approximation, while taking into 
account the quantizer operation. Then, we provide a state 
evolution formalism for analyzing the BP algorithm in the 
large-system limit, and also several theoretical results on the 
impact of finite receiver resolution that can be drawn from 
it. We note that a similar problem has been considered in 
||6l , however based on the replica method from statistical 
physics Q. In (H, linear systems with general separable 



output channels have been considered in the large-system 
limit. Although [8| could include our quantized MIMO case, 
only sparse systems have been considered and no simula- 
tion results have been provided to validate the theoretical 
results. In this work, the "loopy" BP algorithm operating on 
quantized dense linear systems is studied theoretically and 
experimentally. Moreover, our derivation steps for the large- 
system limit are quite straightforward and well justified. The 
main advantage of the BP approach compared to J6) is that 
it is more intuitive and allows to find efficient algorithms 
and analyze their performance and convergence behavior. In 
order to ease calculations, we restrict ourselves to real-valued 
systems. However, the results can be extended to the complex 
case. 

Our paper is organized as follows. Section [TT] describes 
the system model. In Section Hill an approximative BP-like 
detection algorithm operating on quantized data is derived; 
then, we provide a state evolution analysis of the BP algorithm 
and study the effects of quantization in the large-system limit 
in Section [IVJ Finally, in Section [VJ some simulation results 
are presented to numerically validate the theoretical findings. 

II. System Model 

We consider a point-to-point MIMO channel where the 
transmitter employs K antennas and the receiver has N 
antennas. Let the vector x £ M. K comprises the K transmitted 
i.i.d. symbols, each drawn from a certain distribution qo(x) 
with zero mean and variance c x . The unquantized (analog) 
output vector y £ R w is related to the input as 



y 



Hx 



V, 



(1) 



where H £ M. is the channel matrix assumed to be 
perfectly known at the receiver, and r\ refers to Gaussian noise 
vector with covariance R m = E[r]r) T ] = <7qI. 

In a practical system, each receive signal component yi, 
1 < I < N, is quantized by a 6-bit resolution scalar quantizer 
(A/D-converter). Thus, the resulting quantized signals read as 



n = Q(yi), 



(2) 



where Q(-) denotes the quantization operation. For the case 
that we use a uniform symmetric mid-riser type quantizer [9], 
the quantized receive alphabet for each dimension is given by 



k)A; k = !,■■■ ,2 b } =n, (3) 



where A is the quantizer step-size and b the number of 
quantizer bits, which are set the same for all the quantizers. 

With these definitions, the conditional probability distribu- 
tion of the quantized output given an input x reads as 
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Pr\ x (r\x) = Y[p (n\hJx), 



(4) 



1=1 



where h, is the Z-th row of H and 
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with $(a;) represents the cumulative Gaussian distribution 
given by 



1 



exp 



dt. 



Hereby the lower and upper quantization boundaries are 
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,., _A for n >_A (2 b_ 2) 

otherwise, 
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and 

,'n + f for n<f(2 6 -2) 
'-co otherwise. 

III. Approximative BP Detection 

Our goal is to derive a low complexity detector computing 
the conditional mean estimate 



x = E[x\r] 



(7) 



based on the knowledge of qo(x), po(-\-) and H. This problem 
is related to the problem of finding the marginal probabilities 
p Xk \ r (x k \r), for which belief propagation can provide low- 
complexity approximations. To this end, a factor graph repre- 
sentation is needed. 

A. Factor Graph Representation 

In analogy to [4], a factor graph representation of the 
quantized MTMO system is shown in Fig.Q] Each data stream 
Xk is represented by a circle, referred to symbol node, and 
each received quantized signal n corresponds to a square, 
called the signal node. Each edge connecting k and I represents 
the corresponding gain factor hik, if h\% ^ 0. Ignoring the 
cycles in the graph, let us derive the so called "loopy" BP 
algorithm (or sum-product algorithm) from the factor graph 
representation!]] 

'We note that the BP is optimal for cycle free graphs and performs nearly 
optimal in sparse graphs. In the case of dense, large enough, channel matrices, 
it may provide good approximate posteriors as we will see later. 
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Fig. 1 

Factor-graph representation of the quantized MIMO channel. 



B. BP based Detection 

Each iteration t of BP consists in, first sending messages 
from each signal node I to each symbol node k (horizontal 
step), and then vice versa (vertical step). The messages contain 
the extrinsic information of Xk in form of density functions 
computed based on the previously received messages. We 
denote the symbol-to-signal messages by rf k (xk) and the 
signal-to-symbol messages by \\ k (xk): 
Horizontal step: 

A*fcOfc)= / Po{n\hJx) Y[ [Trlk'M&Xk'] (8) 

k'^k 

Vertical step: 

T^Ofe) = otmqo{x k ) Y[ \j k (x k ), (9) 
i'^i 

where aik is a normalization factor so that ir^ 1 (xu) is a valid 
density function. The algorithm is initialized by 

7r ; fe (a; fe ) = q {x k ), l = l,...,N, k = l,...,K. (10) 

C. Approximative BP based Detection Algorithm 

In the case of dense matrices, the complexity of the BP 
algorithm grows enormously with K, since a (K — 1)- 
dimensional integration (or summation) has to be performed 
in the horizontal step. Thus we discuss now an approximation 
scheme, in analogy to HI (see also iflOll ). that may be justified 
in the large-system limit. For this we split the quantity hj x 
as 

hjx = ^ h ik>Xk> + hkXk = wik + hikXk- (11) 

fc'^fe 



Given the distribution p Wik (wik), ® can be rewritten as 

/•oo 

tfk( x k)= / Po{ri\wi k + hi k x k )p Wlk (wi k )dwi k - (12) 



The key remark is that wi k is a weighted sum of the indepen- 
dent random variables x k > distributed according to nj k , {xy). 
Due to the central-limit theorem, it can be regarded as a 
Gaussian random variable in the large-system regime, having 



the mean and the variance 



k'^k 



Hk' > 



Qfc — £ hf k ,v lk ,, 

k'^k 



(13) 
(14) 



respectively, where m\ k , and VA., are the mean and the 
variance of xy according to the distribution n\ k , (xk')- The 
horizontal step can be now written as follows 



tfk( x k) 



Po{n\wik + hikXk) 



-dw[ k 
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n" P ~ /4 ~ hikXk \ _ f r\ ow - p\ k - h lk x k 



Ik 



=Plk(n\l4k + hikXk) 



(15) 



Obviously the approximate BP is a kind of modified parallel 
interference cancellation (PIC) as in the unquantized case 
H), where p\ k represents an estimate of the interference 
component in n for the stream k, and C\ k quantifies its 
MSE. Even if this approximate BP iteration holds in the limit 
of infinite number of antennas, it usually performs well for 
systems of moderate sizes. 

IV. State Evolution Analysis 

State evolution analysis (also known as density evolution 
for the case of sparse matrices) is a powerful tool to study the 
behavior of belief propagation in the large-system limit ifTTl . 
The large-system limit means that we consider the limit when 
K and N go to infinity, while the ratio (3 = K/N is kept fixed. 
Under the conjecture that the presented BP based detector 
would be asymptotically optimal, this analysis would deliver 
useful theoretical results about the MIMO system performance 
under quantization. For the analysis, we assume a random 
channel matrix H, where the entries {hi :k } are i.i.d. with zero 
mean and variance 1/N. The main idea is to approximate the 
messages by Gaussian densities, which holds exactly in the 
large-system limit. In fact, given that hikXk in (TT3T > scales as 
1 / VN, it becomes small as TV becomes large, and as such we 
take the second-order expansion of the messages ( TT3T > as 

Plk( x k\Pik + h ikX k ) «/4(xfc|/4fc) + p\k{xk\l4k) h ikXk+ 



Uk(xMhf k 4 + 0(N-^) 



(16) 

where pjkfaklpjk) an ^ Pik( x k\p*k) denote the first order and 
the second order derivatives of p\ k (xk \p\ k ) with respect to pt . 
Keeping the terms up to the order of TV -1 , and using now the 
following approximation 

1 + ax + hax 2 = e * B +s(6-« i V + Q(x 3 ), (17) 

the horizontal and vertical steps can be represented as 



>Hk( x ) « exp 



'Ik- 1 



it,2 
'Ik 



, and 



(18) 



TT lk (x) oc q (x) exp 



qt,2 

'l'k 



I'^l 



(19) 



respectively, where we introduced the definitions 



Ulk Piking) k M 

For the following, we assume a random matrix H with i.i.d. 
entries and we consider the large-system limit. The state evo- 
lution analysis aims to study the dynamics of the BP detector, 
i.e. to track the evolution of the densities parameters, namely 
p\ k , C\ k from ( fT3] l and (fl4l i. respectively, and additionally 

z ik — ®v k > $ik — @vk ' G\ k = £ a*/ k (21) 

v^i I'^ti i'jti 

over the iterations. Strictly speaking, these variables are ran- 
dom. Relying on a heuristic assumption that the incoming 
messages (6\ k , S* fc ) to each symbol node remain independent 
from iteration to iteration, when conditioned on a given 
transmit symbol Xk and channel H, and by the central 



and 



(20) 



limit theorem, we conclude that p,j k , 



(from ( fT3l l and 



(O), z\ k , Ff k and G\ k become asymptotically Gaussian. In 



particular C\ k , F\ k and G\ k are sums of terms of order 1/N, 
thus they admit asymptotically zero variance. In other words, 
they become deterministic in the large-system limit at given 
iteration provided that H has i.i.d. entries, i.e. independent 
of the given realization of H, the received vector r and the 
indexes I and k; that is 

hfk'Vi, -> Yl K^ik'V^lH^k}^, (22) 

(23) 
(24) 



k'^k 



k^J2 



^ rWh [9Y k \H,x k ]=F\ 



£"i'fc ~* ["f'fcl-H', Xk}— G\ 

where — > symbolizes the convergence to the asymptotic limit 
Therefore, we can write the messages rf k (x) as 



Tf lk (x) -> aq (x) exp 



z\ k x-\{F t -G t )x 2 



(25) 



We will return to the calculation of these parameters later on. 
Now let us consider the joint distribution p WlhllJ t (w, p) of the 
interference term in antenna I for symbol k, wi k from ( fTTI ). 
and its estimate p\ k in ( fT3l , given H and x k . Again by the 
central limit argument, wik and p5 k are asymptotically jointly 
Gaussian 



H,xi 



with the covariance matrix R w , having the entries 

E[wf k \H,x k ] =E[( ]T h lk' x k')' 

k'jik 

= ^2hl,E[x 2 k/ \H] 

k'^k 



(26) 



/3 / x k >q (xk>)dx k >=l3c x 



E[pf k 2 \H,x k ] =E[(£ %,mk,) 2 |H,a*] 

k'^k 



E[w; fe /i* fc |ff, ifc] = ^fc//ijjfe»E[s k »mf k /|ir] 

= £ h 2 lk ,V[xh>™\ k ,\H] (28) 



/3Ez . 



Xk' / ^7T*fe' (x)dxC 



In the steps above, we used the fact that the expectations with 
respect to xy and z\ k , become asymptotically independent of 
the indexes I and k'. Thus, we obtain the bivariate Gaussian 
distribution 



— M)" 

"f72 — 



2 H (<v,, ; 



Pw lh ,n\ k {w,p) = 



(29) 



Next, the MSE parameter C* in ( l22l can be expressed as 

C 4 = ^ E [h 2 lk , f x^\ k , (x)dx ~ rr$\H, x k ] 



k'^k 



--J2hf k ,E[ x 2 n t lk ,(x)dx H,x k \-J2hl>n™j> k 2 \H] 

k'^k J k'^k 



x 2 7r t lk ,(x)dx -(301=13(4 -c* m ). (30) 



Note that the densities p\ k {-\-) defined for the horizontal step 
in ([T5l l become independent of I and k, and read as 



p*(r\n) = $ 



r low _ ^ 



(31) 

Afterwards, we compute F t from ( |23l . For that we need the 
joint distribution p r ti tj Xk (r,p t [k \x k ) f° r f" xec l channel H 



Pr 4 i\ k \x h ( r ,^lk\ x k) = J Pr,^ lk ,w lk \x k ( r p\ k ,Wl k \x k )dwi k 

=Pr l \ Xk {ri\x k ) J p rh ^ ikiWlklxk (ri,f4 k ,wi k \x k )dwi k 
=P ri \x k (ri\x k ) p (ri\wi k + h lk x k )p Wlk ^t^w lk , ri k )dwi k 



-+Pn\x k {n\%k)p t (n\t4k + h ikx k ) 



t,2 
23c*,, 



(32) 



where ri is the vector containing the elements of r excluding 
77 and the density function /?* ( • | • ) is obtained using d29l by 



performing the integration as 



(27) p*(r|^) = $ 



_d> 



r low _ ^ 



'/3(c x -^)+a ( 2 





(33) 

From ( 1231 . (f20b and d32b we identify after some straightfor- 
ward steps the asymptotic non-vanishing term for the param- 
eter F t 



F 



p\r\p) 



pHr\p) 



2/3c 



-.dp, (34) 



where we dropped the indexes I and k due to asymptotic 
independence in the large-system limit. Similarly, we obtain 
the parameter G* in d25T l as follows 



G 



E 



p\r\p) 



j?{r\p) e ^ 



P*{r\n) 



= d//. (35) 



We turn now to determine the distribution of z\ k defined 
in ( ETl i. which, as mentioned before, follows a Gaussian 
distribution conditioned on x k . Again by (f20b and d32b . we 
show that its mean is 



^r\x k ,H[Zlk\ x k) H] 



,P t ( r \f J ') 



hv 



2/3c 



I'k- 



zdp 



E p^f-t 



P*(HaO e ^ 



-.dp ■ x k =E l ■ x k , 



(36) 



P*(r|M)V2^4 
and its variance is 

E r |i fc ,ff[^fc 2 kfe I -H'] - ( E r|x fc ,H'[^fekfe ) - H ']) 2 

^E E H^^>^] ift ' (3?) 

I'^l 

Note that z* fc « ^f=i fy'fc = 4 and thus 7r* fc (x fe ) « tt* (x k ) 
for all 2 because each 9i k has a vanishingly small effect on 
the sum. In summary, we get the conditional density 



Pa(4\ x k,H) = 



1 



r,t pt T -,2 
(z,-E m k ) 



y/2nF t 



(38) 



Since the BP iteration is initialized by ni k (x k ) = qo(x k ), 
we have from $T7\ , (l28l and (l30l the initial parameters c" t = 
: and c? = c x , it can be shown by mathematical 



,,0 



induction that for all t 



-~x,rm ^x 



c x , p\r\p) = p\r\p), 



& = F*, G* = 0, 74 +1 (z fe ) oc 9o (zfc) -pUzk^k)- 



(39) 



Finally, we conclude that the performance of the large-system 
regime is fully described by the sequential application of two 
updates for the parameters F t and c\ n (cf. (l34l and d27l >) 



F 



E 



ren 



p'(r| M ) 2 



[/ xpc(z|a;)go(^)dx 



=d/i, 



(40) 



-dz, 



Jp^(z|x)g (a;)dx 

with the initial state = 0. Interestingly, the stationary 
conditions of the state evolution equations coincides with the 
fixed point equation found in |7| with the replica method. 

V. Numerical Results 

Let us consider BPSK transmission, i.e. Xk G {— 1, +1}. 
We have c x = 1 and from d40b we get 



J+ 1 — 



1 



2V2ttF 
1 



2f* — e 2F* J tanh(z)dz 



e 2ft tanh(z)dz, 



(41) 



where the second line follows from the antisymmetric property 
of the tanh(z) function. The equations d40b of the state 
evolution as well as the typical trajectory of the BP iteration 
are shown in Fig. [2] for (3 = 1.8, Co = 0.1 and 6 = 4. 
The quantizer step size A has been chosen to minimize the 
distortion under Gaussian input. We distinguish different fixed 
points. Clearly, the performance of BP based detection algo- 
rithm is characterized by the poor solution as shown in Fig. [2] 
Besides, from the distribution of given Xk, p t G {zk\xk) in 
(l38T l. which directly affects the decision of Xk through iTk(xk), 
it immediately follows that the bit error probability (BER) after 
performing t iterations is given by 

P l b = Pr(^ < 0) = Q(VF*), (42) 

where Q(z) = f°° —^c^^dt is the error function. The 
analytical BER performance at the fixed point for j3 = 1 (N = 
K) is shown in Fig. [3] as function of the SNR= 10 log 10 \ 
for different number of bits. The experimental BERs, carried 
out for a 20 x 20 system with b e {1,4} and using 10 BP 
iterations, are also shown for comparison. Obviously, there is 
a good match between the theoretical results and the Monte 
Carlo results for 6 = 1, while there is a gap for b — A due 
to the insufficiently high number of antennas. Nevertheless, 
the analytical curve for b = 4 still predicts correctly that 
the performance loss with four bits compared to the ideal 
case (b = oo) becomes negligible. We note that d4Qb for 
the noiseless case (do = 0) always admits a possible perfect 
detection solution, i.e. = 1 and F°° = oo [6|. However, 
since the fixed point solution is not unique and BP usually 
converges to the worst solution, the BER behavior over the 
SNR might exhibit an error-floor as shown in Fig. [3] for the 
case 6 = 1. The minimum number of bits needed to ensure 
perfect detection in the noiseless case, i.e. where the recursion 
(l40l evolves to the unique fixed point solution at = 1, 
is depicted in Fig. |4] as function of the load factor f3. We 
observe that for a system load f$ < 1 even 2-bit ADCs might 
be sufficient for symmetrical systems. 




Fig. 2 

State evolution chart, /3 = 1.8, ctq =0.1 and b = 4. 
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Fig. 3 

1 AND DIFFERENT BIT 
RESOLUTIONS. THE SIMULATIVE RESULTS WERE OBTAINED FROM THE 
APPROXIMATIVE BP DETECTION APPLIED ON A 20 X 20 SYSTEM. 



VI. Conclusion 

We studied a low complexity detection algorithm based on 
belief propagation for quantized MIMO systems. Additionally, 
a state evolution formalism has been presented to analyze the 
performance of the BP based detector when operating on quan- 
tized data in the large-system limit of fading channels. A set of 
simulation results was provided, which, in agreement with the 
analytical results, shows that the BP approach achieves good 
performance even with low resolution ADCs and moderate 
number of antennas. 
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